Feedback control of the squeezing of the 
fluorescence light 

Alberto Barchielli and Matteo Gregoratti 
Politecnico di Milano, Dipartimento di Matematica, 
Piazza Leonardo da Vinci 32, 20133 Milano, Italy 
and INFN (Istituto Nazionale di Fisica Nucleare) 

October 12, 2009 

Abstract 

Among the formulations of the theory of quantum measurements in 
continuous time, quantum trajectory theory is very suitable for the intro- 
duction of measurement based feedback and closed loop control of quan- 
tum systems. In this paper we present such a construction in the concrete 
case of a two-level atom stimulated by a coherent, monochromatic laser. 
In particular, we show how fast feedback a la Wiseman and Milburn can 
be introduced in the formulation of the theory. Then, the spectrum of 
the free fluorescence light is studied and typical quantum phenomena, 
squeezing and sub-natural line-narrowing, are presented. 

1 Introduction 

The theory of continuous measurements and filtering for quantum systems [6J [1] 
has opened the possibility of a consistent introduction of measurement based 
feedback in continuous time and, so, of closed loop control of quantum systems. 
The formulation of continuous measurement theory adopted in the paper is the 
one based on stochastic differential equations, commonly known in quantum 
optics as quantum trajectory theory [7j. 

Photo-detection theory in continuous time has been widely developed inside 
quantum trajectories [3] and applied, in particular, to the fluorescence light 
emitted by a single trapped two-level atom stimulated by a coherent monochro- 
matic laser. As well as various feedback schemes on the atom evolution, based 
on the outcoming photocurrent, have been proposed [9]. However the introduc- 
tion and the analysis of feedback have been mainly focused on the control of 
the atom. The typical aim was to drive the atom to a preassigned asymptotic 
state. 

Here, we are interested not only in the atom, but also, and mainly, in the 
emitted light. Of course, a detection of the emitted light can be seen as a direct 
observation of the emitted light as well as an indirect observation of the emitting 
atom mediated by the fluorescence light. Our aim is to employ control and 
feedback processes to enhance the squeezing properties of the fluorescence light. 
The squeezing of the fluorescence light can be checked by homodyne detection 
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and spectral analysis of the output current. We consider the mathematical 
description of photo-detection based on quantum trajectories, as it is suitable 
both to consistently compute the homodyne spectrum of fluorescence light, and 
to introduce feedback and control in the mathematical formulation. We study 
how the squeezing depends on the various control parameters, how feedback 
mechanisms can be successfully introduced and used to enhance the squeezing 
[H [3] . We consider only the Markovian feedback scheme a la Wiseman-Milburn 
[S], which is simple, but flexible enough to give a physically interesting model 
of closed loop control. We also show that, under certain conditions, feedback 
can produce an effect of line narrowing in the homodyne spectrum, similarly to 
what happens to an atom stimulated by squeezed light. 

Apart from the new results on enhancing squeezing by feedback and on line 
narrowing, there is a conceptual interest in using quantum trajectory theory 
in deriving the homodyne spectrum. The traditional approach to homodyne or 
heterodyne spectrum is to define, by some analogy with the classical case, a suit- 
able quantum correlation function for the outgoing electromagnetic field. Then, 
the spectrum is defined to be the Fourier transform of such a correlation func- 
tion. Finally, the quantum fields are eliminated in favour of the atomic variables 
and the "quantum regression theorem" is used to get the final result. On the 
other side, quantum trajectory theory is based on quantum measurement theory 
and gives the output of the measurement in continuous time together with its 
distribution, consistently with principles of quantum mechanics. As the output 
is a stochastic process, its spectrum can be rigorously defined by probability 
theory by using its distribution. Moreover, thanks to the consistency of quan- 
tum continuous measurement theory, even when the quantum fields have been 
already eliminated (so that quantum regression theorem is implicitly already 
contained), the output can be seen not only as an observation of the system, 
but also as the result of an observation of the output field which mediates the 
measurement. To succeed in describing a quantum effect such as squeezing of 
the fluorescence light by using quantum trajectory theory is to show that such 
a theory, in spite of its "classical flavour" , is fully quantum mechanical. 

We consider a trapped two-level atom with Hilbert space C 2 . Denoting by 
a the vector (a x , <7 y , a z ) of the Pauli matrices, let the free Hamiltonian of the 
atom be loqg z /2, with resonance frequency luq > 0, and let the lowering and 
rising operators be cr_ and er + . Let the eigenprojectors of a z be denoted by P+ 
and P- and, for every angle 4>, let us introduce the unitary selfadjoint operator 

= e 1 ^ er_ + e^ 1 ^ a + — cos 4>cr x + sin </> a y . 

2 The dynamics 

We start by giving the theory in the case of no feedback and then we show how 
to introduce it, according to the scheme of [S]. 

2.1 No feedback 

We admit an open Markovian evolution for the atom, subjected to 'dephasing' 
effects and to interactions both with a thermal bath and with the electromag- 
netic field, via absorption and emission of photons. The atom is stimulated by a 
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coherent monochromatic laser. We denote by 7 > the natural linewidth of the 
atom, by Q > the Rabi frequency by u> > the frequency of the stimulating 
laser and by Aw = ujq — u> the detuning. Other parameters are the intensities of 
the dephasing and thermal effects, k^ > and n > 0. The state n t is governed 
by a Master equation, which turns out to be time-homogeneous in the rotating 
frame (which we adopt here and in the rest of the paper): 



where 



dt 



m = Cm , 



Alu 



— a z + -a x ,p 



+ 7^d (<?zP<?z - P) 



( <7+p<7- - ~ {P- , p}j + 7(n + 1) (a- pa + - ~ {P + , p} 



Note that the intensity of the laser enters only in the effective Hamiltonian 
through the Rabi frequency f2. 

The measuring apparatus is made by two homodyne detectors. Part of 
the emitted light reaches the detectors and part is lost in the free space. The 
fraction of light detected by one of the detectors depends on its efficiency, on the 
spanned solid angle and can eventually be enhanced by using a focussing mirror. 
So, the fluorescence light is divided in three parts according to the direction 
of propagation: we call side channel k (k — 1, 2) the directions reaching the 
detector fc, and forward channel those of the lost light. The stimulating laser 
is well collimated in such a way that it does not hit the detectors; so, we can 
say that it acts in the forward channel. We denote the effective fractions of 
light emitted in the forward and in the two side channels by |ao| 2 , |ai| 2 , |ofa \ 2 , 
respectively; obviously, |ao| 2 + |cti| 2 + \a2\ 2 = 1. For k = 1,2, we can also 
interpret |afc| 2 as the total efficiency of the detector k. Moreover, the initial 
phase of the local oscillator in each detector is denoted by -&k and it is included 
in the parameter at € C by setting = arg ctk ■ To change means to change 
the measuring apparatus. 

We can condition the evolution of the atom on the continuous monitoring 
of the photocurrents. The two homodyne photocurrents 1\ and I2 and the con- 
ditional state of the atom p t (a posteriori state) are stochastic processes whose 
distributions depend on the initial state of the atom and on the parameters in- 
troduced so far. In particular the atom has still a Markovian evolution, even if 
stochastic. Let us introduce the latter by means of the linear stochastic Master- 
equation 

2 

da t =£o- t dt + y /j^K[a j a-]<TtdW j {t), (1) 
j'=i 

where, for every matrix a, the superoperator 72. [a] is 

TZ[a]p = a p + pa*, 

and where W\ and W2 are two independent standard Wiener processes in some 
reference probability space(f2, J 7 , Q). In all these equations the initial condition 
is a statistical operator Qo = rjQ = o~o and the solution <j t of U} is a positive 
operator valued stochastic process. 
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Then, the homodyne photocurrents are the (generalised) stochastic processes 
Ij(t) = Wj(t), while the a posteriori state is 

p t :=—, with p t = Tr{cr t }. 
Pt 

The important point is that the physical distribution of the homodyne currents 
and the a posteriori states in the time interval [0, T] is given by the new prob- 
ability Fxidoj) = pt (w)Q(do;): the quantity pr = Tr{or} is the density of the 
physical probability with respect to the reference probability. By the fact that 
Ptit>0, turns out to be a Q- martingale, we have that the physical probabilities 
do not depend on the final time T. 

Of course, we could switch on the detectors, but decide to ignore the outputs. 
This should not modify the evolution of the atom and, indeed, for every t < T, 

T) t = t a t (cj)Q(duj) = f pt(u)F T (du). 

Thanks to the linearity of (TTJ) , we can introduce the stochastic evolution map 
(or propagator) A(t,s), satisfying 

2 

A(t + dt, t) - Id = Cdt + V7$I K[aj<r-]dWj(t), 

3=1 

so that o~t — A(t, 0)<To. Here, Id is the identity map on the space of 2 x 2 
complex matrices. 

2.2 Introduction of the feedback 

We introduce a feedback scheme based on 1%, in order to modify the properties 
of the fluorescence light in channel 2. We check such properties by analysing the 
properties of I2, but, of course, we could remove the second homodyne detector 
and employ the light emitted in channel 2 for other purposes. Let us call it 
free light. Our scheme is summarised by the following picture, where "h. det." 
means "homodyne detector" : 




laser 



Assuming instantaneous feedback, the amplitude of the stimulating laser is mod- 
ified by adding a term proportional to I\ , with the same frequency u and with 
initial phase possibly different from that of the original laser. Let this phase 
difference be <p € [0, 7r). As the laser intensity appears only in the Hamilto- 
nian part of the evolution (TTJ) , the effect of the feedback is to give rise to a new 
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Hamiltonian term I\{t) M with M = c^/ja^, where eel controls the intensity 
of the feedback. Let us deduce the modified evolution equation. By defining 
the map Mp — — i[M, p], the contribution of the feedback to the propagator in 
an infinitesimal interval turns out to be 

e MdWl(t) =ld + MdW 1 (t) + ^M 2 dt. 

Taking into account that the feedback must act after the signal is produced, the 
new infinitesimal propagator is 

Ai(t + dt,t) =e MdW ^ oA(t + dt,t). 

By the Ito rules we obtain 

A f (t + dt,t) - Id = Cfdt 

+ VT^K o- - iccTy] dWi(i) + V^n[a 2 a-] dW 2 {t), 



where 



dp = 



<?z + tt °x , P 



2 2 



+ 7 fc d (pz p<t z - p) 



-11 [ (j+pa-- , p} 



+ 7(n+ 1 - |ai| 2 ) (o-pcr+ - \{ p + , P} 
+ 7(01 cr_ - \ca v ) p (a i a + + ica v ) 
- |{(K| 2 - 2c| ai | sin(l?i - if)) P+ + c 2 , p} , 

Aw c = Auj + c-f |ai| cos(i?i —ip)&R. 

Thus, thanks to the features of the particular feedback introduced, we get an- 
other linear stochastic Master equation 

dcr t = Cfcrt dt + VY^N a- - ica v }a t dWi(i) + ^U[a 2 o_\o t dW 2 (t). (2) 

By Girsanov theorem, a fundamental result of stochastic calculus, it is possi- 
ble to prove that under the physical probability the output homodyne currents 
can be written as 

Ij{t)dt = dWj{t) = dBj(t) + Vj(t)dt, 

where £?i(t), B 2 {t) (t € [0, T]) are two independent standard Wiener processes 
under Ft and 

vi (t) = 2^ReTr{{a 1 a--ica ip )p t }, v 2 (t) = 2^Re (a 2 Tr {a_p t }) . 
Then the mean function of I 2 is 

E PT [I 2 (t)]=v 2 {t). 
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Moreover, by using techniques based on characteristic functionals, explicit ex- 
pressions for the higher moments of the output have been obtained [3]. In 
particular, the auto-correlation of I2 is 

Ep T [I 2 (t) h{s)] = <S(t - s) + 7 Tr [n[a 2 cT-} 

o e |4 - s|£ ^ o n[a 2 <T-} o e min ^ c ' g }- (3) 

The expression above is similar to what is obtained in traditional approaches 
through the quantum regression theorem; the general theory of continuous mea- 
surements guarantees that it is a positive definite function of s and t as it must 
be for an auto-correlation function. 



3 The spectrum of the free light in channel 2 



Let us consider now the spectrum of the light in channel 2, the fluorescence light 
detected but not involved in the feedback loop. The homodyne current 12(f) 
is a stochastic process and it can be proved that it is asymptotically station- 
ary. Then, the spectrum S% (fj,) is the Fourier transform of the autocorrelation 
function for long times. This classical definition can be recast in the following 
way: 

S 2 (fi) = Sel(fl) + SinelC/u), 



SdO) = lim h 

I — *+oo 1 



e^ s I 2 (s) ds 



SinelO) = lim 7f 



e ifls I 2 (s) ds 



I e^ s I 2 (s)ds 
Jo 



>0, (4) 



where the spectrum is decomposed in the elastic and inelastic parts. By explicit 
computations on can see that S c i(/i) is proportional to a Dirac delta centred in 
zero (as we are working in the rotating frame, /1 = corresponds to resonance 
with the frequency of the stimulating laser). Instead, S mc i(^i) is the limit of the 
normalised variance of the Fourier transform of the photocurrent I2 ■ An explicit 
expression can be obtained from the expression of the second moments ([3]) . The 
asymptotic behaviour of the atomic a priori state r\t = e Cft Qo ensures that the 
limit defining the spectrum exists and that it is independent of the initial state 
go of the atom. We get 



SinelO) = 1 + 2j\a 2 \ 2 S ■ 



A 



A 2 



(5) 



where s = (cos d 2 , sin -d 2 , 0) , the matrix A has matrix elements 
ai3 = a 3 i = 0, a 23 = -a 32 = fi, 
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Figure 1: S ine i(fi) with 7 = 1 and: (1) Acu = 0, c = 0, ft = 0.2976, i? 2 = -tt/2; 
(2) Aw = 1.8195, c = 0, ft = 1.7988, i? 2 = -0.1438; (3) Alo = 0, c = 0.0896, 
17 = 0.2698, 0! = tt/2, i9 2 = -tt/ 2 > <ys = 0; (4) Aw = 1.6920, c = 0.1326, 
ft = 1.9276, i?i = 2.8168, d 2 = -0.0851, <p = 1.2460. 



an = 9 ^ ^ d 2c l ai I cos ^i sm V + 2c 2 sin 2 
ai2 = Auo c — j(c\ai\ cos(0i c 2 sin2<£>J , 

a 2 i = — Auj c — 7^c|ai| cos(??i + <p) + c 2 sin 2^^ , 
a 22 = 7^- + n + 2kd — 2c|«i| sin$i coscp + 2c 2 cos 2 <pj , 

a 33 =7(1 + 2n- 2c|ai|sin(i?i - ip) + 2c 2 ), 
f = Tr I (e" 52 cr_ ?7 oq + e~ I,5fc ry cq cr + - Tt{<7$ 2 r) cq } r) cq ) a\. 

and ?7 q = lim e ft Qo- 

In the case |a 2 | = we get Si ne i(fi) — 1, which is the spectrum of a pure 
white noise. Indeed in this case no fluorescence light reaches the detector and 
we see only the spectrum of the fluctuations of the local oscillator, interpreted 
as shot noise. 



3.1 Squeezing 

In general it is possible to give a description of the output J 2 in terms of a 
measurement of some Bose quantum fields, to show that the value of Si n el(/-*j $2) 
is the variance of a quadrature of the field in channel 2, the value of Si ne \(p; $2 + 
7r/2) is the variance of the conjugate quadrature 2 . Then, the Heisenberg 
uncertainty relations imply that 

Si„ei(£i; 2 ) S iaB i(ji; 2 + 7t/2) > 1. 
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If Si nc \(p; $2) < 1 for some /x and 62, the field is said to be squeezed. 

Some typical spectra, showing well pronounced squeezing and produced by 
suitable choices of the control parameters, are given in Figure [TJ The position 
of the minimum can be tuned by a suitable choice of the control parameters; 
in Figure [T] the minima are in fj, = and in fi = ±2.5, with and without 
feedback. The lines (1) and (2) are without feedback and the lines (3) and 
(4) are with feedback. Here and in all the other graphical examples we take 
|ai| 2 = \a 2 \ 2 = 0.45, n = and k d = 0. 

Note that the Heisenberg uncertainty relations give rise to peaks in the 
complementary field quadratures, as shown in Figured 




-6 -4 -2 2 4 6 



n 

Figure 2: S- m ei(fi) with all the parameters as in the previous figure, but $2 
tf 2 +7r/2. 



3.2 Line-narrowing 

By feedback control another quantum effect can be produced, the line-narrowing. 
After the first observation of squeezing, Gardiner predicted that stimulating a 
two-level atom with squeezed light would inhibit the phase decay of the atom. 
The squeezed light would break the equality between the transverse decay rates 
for the two quadratures of the atom and one decay rate could be made arbitrar- 
ily small, producing an observable narrow line in the spectrum of the atomic 
fluorescence light. This was seen as a "direct effect of squeezing" and thus as a 
measure of the squeezing of the incident light. Nevertheless, Wiseman showed 
that this atomic line-narrowing is not only characteristic of squeezed light, but 
it can also be produced by immersing a two-level atom in 'in-loop squeezed' 
light. The difference is that in the Gardiner case the other decay rate becomes 
larger, while in the Wiseman case it is left unchanged. 

Now we can show that the same atomic line-narrowing can be obtained 
stimulating a two- level atom even with non-squeezed light, that is with a coher- 
ent monochromatic laser in presence of a (Wiseman-Milburn) feedback scheme 
based on the (homodyne) detection of the fluorescence light. This effect is ob- 



8 



- 4c =0.45 
-C=0 




Figure 3: S ine i(fi) for c = |a x | /2 or c = and fi = 2, p x = p 2 = 0.45, Alo = 0, 
v = 1,^=^,^ = 0, 7=1. 



tained within the model described up to now, in a different region of the control 
parameters. 
Let us take 

<p = | , 0i = 7T, ^2 = 0, Acj = 0, n = 0, fc d = 0. 

Then, in the case of no feedback c = one obtains a spectrum which is a 
Lorentzian peak with width 7: 



2|a 2 | 2 f7 2 - 



5inel(M) = 1 + 



7 2 + 2ft 2 /i 2 +7 2 /4 



Instead, with the optimal choice of the feedback, c — \a\ \ /2, we get a Lorentzian 



of width I 1 — I ax I ) 7 



M 2 (4^7 2 + 2» 2 ) (i-M 2 ) 7 2 

*Sinel(Mj — 1 H 7 : ; ~~o 



This is the sub- natural line- narrowing effect, shown in Figure [3] 



4 Conclusions 



The new results presented either are of general conceptual relevance, either 
clarify the behaviour of fluorescence light and atom under coherent driving and 
feedback. 

First of all we have shown how to introduce the spectrum in quantum trajec- 
tory theory through its classical definition ^ , without ad hoc quantum defini- 
tions, but in agreement with the axiomatic structure of quantum measurement 
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theory and with probability theory We have verified that this agreement with 
quantum measurement theory ensures consistency with the existence of the 
traced out quantum field, so that Heisenberg uncertainty relations hold for such 
spectra and the squeezing of the output field can be analysed. 

We have shown how feedback (at least in the form a la Wiseman-Milburn) 
modifies the spectrum of the free fluorescence light and can enhance its squeez- 
ing. Many physical effects are introduced all together: detuning, thermal and 
dephasing effects, not perfect detection efficiency, control. The final formula for 
the spectrum (O is given with all the parameters introduced by these effects. 

Finally, we have shown that feedback can produce line-narrowing in the free 
fluorescence light, even if the atom is not illuminated by squeezed light, but 
only by coherent light. 
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